Normal Forms of Hopf— Zero singularity 



Majid Gazor*, Fahimeh MokhtariO 



Department of Mathematical Sciences, Isfahan University of Technology 
Isfahan 84156-83111, Iran 
December 11, 2012 



Abstract 



The Lie algebra generated by Hopf-Zero classical normal forms is decomposed into 
two versal Lie subalgebras. Some dynamical properties associated with each subalge- 
bra are described; one is the set of all volume-preserving conservative systems while the 
other is the maximal Lie algebra of nonconservative systems. This introduces a unique 
conservative-nonconservative decomposition for the normal form systems. There exists a 
Lie-subalgebra that is Lie-isomorphic to a large family of vector fields with Bogdanov- 
Takens singularity. This gives rise to a conclusion that the local dynamics of Hopf-Zero 
singularities is well-understood by the study of Bogdanov-Takens singularities. Despite 
this, the normal form computation of each singularity is independent of the other. Thus, 
by assuming a quadratic non-zero condition, complete results on the simplest Hopf-Zero 
normal forms are obtained in terms of the conservative-nonconservative decomposition. 
Some practical formulas are derived and are implemented using Maple. The method has 
been applied on the Rossler and Kuramoto-Sivashinsky equations to demonstrate the 
applicability of our results. 

Keywords: Normal form; Hopf-Zero singularity; s[2-representation; Conservative and 
nonconservative systems. 
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1 Introduction 

In this paper we are concerned with computing the simplest normal form of the system 



where h.o.t. denotes nonlinear terms (higher order terms) with respect to {x,y,z) G M?. The 
system (11.11) can be transformed into the first level (classical) normal form 
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X := h.o.t., y := z + h.o.t. 



z := —y + h.o.t. 
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where also [H Equation 2.15]. All the existing results on the simplest normal 

forms of this singularity have only dealt with the cases where 020^10 7^ 0, see [H [3 [9l |26]. The 
main reason for this is that terms associated with 020, foio? Cio? and aoi all have the same 
grades in the usual graded Lie algebraic structures; i.e., the gradings may not distinguish the 
associated four vector fields. Therefore, one needs to deal with all four vector fields in the 
computations which is a tremendously difficult job. A novelty in the results of Algaba et. al. 
P was to notice that the first two non-zero terms (associated with 020 and 610) play the key 
role in the computations and thus, generic conditions with respect to 020 and foio are assumed. 
Thus, they obtained the simplest normal forms and the simplest orbital normal forms. Chen et. 
al. [Bl |9] approached this family of systems using an essentially different method and provided 
an independent solid proof for their uniqueness. Yu and Yuan |26] made an efficient computer 
program to compute the simplest normal form of this family of systems. We do not address 
the convergence of normal forms in this paper. Therefore, all our claims with regards to the 
dynamics are limited to the formal flows of the formal systems. 

A noble goal of normal form computations is for its possible applications in bifurcation 
and local stability analysis. We consider the cases that have not already been considered in 
the literature. Further, our results have other advantages such as more simple computations, 
dynamically meaningful decomposition of vector fields, revealing some possible symmetries, 
possible application in developing computer ODE-solvers, classification and introduction of 
important families of vector fields, descriptive study of the Lie subalgebras, a comparison 
between dynamics of different singularities, and possible applications for more complicated 
singularities. 

A motivation is to use the term associated with aoi as the main player in calculations; this 
is to complete the existing results on this problem with a quadratic non-zero term. (Throughout 
this paper the only assumption is aoi 7^ 0.) The associated grading in our approach distinguishes 
this vector field (grade 1) from the other three (grade 2). This is achieved through a s[2- 
representation for the classical normal forms. This greatly simplifies the computations; we 
obtain complete results on the simplest normal form of these systems without any extra generic 
conditions. Our approach can be applied on many well known models such as Roseler and 
Kuramoto-Sivashinsky equations. 

Dynamically meaningful decomposition of a system is a substantial asset in applying many 
techniques in the study of nonlinear dynamics. A powerful tool for such kind of decompositions 
is s(2-representation of vector fields. This has been mainly applied on nilpotent singularities. 
For example, this decomposes a two dimensional Bogdanov-Takens normal form into sum of 
an Eulerian and a Hamiltonian vector field. This kind of decompositions is greatly appreciated 
in the global analysis of nonlinear dynamical systems such as Melnikov theory, perturbation 
theory for normally hyperbolic invariant manifolds, and Kolmogorov-Arnold-Moser (KAM) 
theory, see p5l Chapter 33]. Furthermore, the dynamics analysis of each decomposed vector 
field helps in a better understanding of the whole vector field as if each component competes 
and interacts with the others in order to determine the whole dynamics. For instance, this has 
been used in the study of reaction-diffusion equations, where reaction and diffusion terms do 
the interactions and competitions to describe the dynamics. 
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There are at least two more advantages for such decompositions stated by [6l [23]. Differ- 
ent sphtting procedures of differential equations can be used in designing geometric integrals 
and decompositions can be used in an ODE-solver by knowing an integral for a piece of the 
vector field. This can also have qualitative and quantitative purposes, see [6]. Barrio and 
Palacian propose an algorithm working with Lie transformations of ODEs that has been split 
into conservative and dissipative components. Their algorithm is stated as more efficient than 
the standard method, since they work with scalar expressions of the conservative part rather 
than its associated vector field terms. 

The computational burden has been the main obstacles of most works in the literature. A 
systematic approach is required for such computations; a Lie-graded structure is an important 
instrument where grade-homogenous parts of vector fields are simplified inductively. In this 
process a basis for the space of grade-homogenous vector fields are used. The mostly used 
choice for the basis has been the monomial vector fields oblivious of their dynamics. This 
usually involves large matrices and it is hard to find the patterns of computations. The con- 
servative and nonconservative polynomial vector fields play this role in this paper instead of 
the monomials. Then, computations yield more solid patterns compared to when monomial 
vector fields are used. Indeed both transformation generators and normal forms are presented 
through conservative and nonconservative grade- homogeneous vector fields. This is a new fea- 
ture of this paper that distinguishes our results from the existing results on the simplest normal 
forms of this singularity. This is accomplished via a 5[2-representation for the classical normal 
form vector fields. This may be surprising to see the theory of s[2-representation is applied to 
a non-nilpotent system. Another novelty of our results is to use sl2-style for the second level 
normal forms; s(2-styles have only been used for the first level normal forms in the existing 
literature. 

Any Lie subalgebra structure generates a group of transformations and the group presents 
a classification within the Lie algebra. These usually introduce important families of vector 
fields. We denote =Sf for the Lie algebra generated by the first level Hopf-Zero normal forms. 
Two transversal Lie subalgebras for ^ are presented. These represent the quasi-Eulerian 
vector fields and volume-preserving vector fields with a first integral. Furthermore, there exists 
a Lie subalgebra from that is Lie-isomorphic to a Lie subalgebra J^b from the Lie algebra 
generated by all two dimensional vector fields with Bogdanov-Takens singularity. This provides 
an explanation for why the s[2-representation works fine for non-nilpotent singularities. This 
also suggests that it may be applicable to some other non-nilpotent singularities. Further, this 
can have an interpretation in terms of their dynamics. 

We conclude that the local dynamics of a Hopf-Zero normal form is well-understood by 
the study of Bogdanov-Takens singularities. The basis of our claim is as follows. There is a 
common practice to ignore the phase component of the Hopf-Zero normal form and obtain a 
planar reduced system. Then, study the dynamics of the planar reduced system and extract 
the full three dimensional dynamics from the planar reduced system. Given this, we prove 
that for any planar reduced system from a Hopf-Zero normal form, there exists a Bogdanov- 
Takens singularity from =Sfb such that the flow of the planar reduced system is embedded into 
the flow of the Bogdanov-Takens singularity, see Theorem 12.71 It is interesting to note that 
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the complexity of Bogdanov-Takens singularities are expected to be more rich than that of 
the planar reduced systems obtained from Hopf-Zero. This is because of two reasons. First, 
the embedding between the two flows assigns the square of the amplitude variable from the 
reduced systems into a state variable from Bogdanov-Takens systems that is not necessarily 
nonnegative. This excludes certain dynamical complexities. Second, a complement space to J^b 
is expected to absorb most of the dynamics of Bogdanov-Takens singularity. This is because 
all nonzero nonlinear terms from =Sfb appearing in a Bogdanov-Takens singularity are simplified 
in the s[2-classical normal forms. This contributes to the exclusion of most complex dynamical 
behaviors. 

The rest of this paper is organized as follows. The conservative-nonconservative decompo- 
sition for normal forms are introduced in Section |2l This is achieved by a s[2-representation 
for =Sf and presenting two transversal Lie subalgebras. Some properties associated with each 
family are described. The second level normal form is computed in Section [31 In Section H] we 
obtain the simplest normal forms. The procedure is divided into three cases and accordingly 
their simplest normal forms is computed in three subsections. Some practical formulas are 
derived and applied on Rossler and Kuramoto-Sivashinsky equations in Section [51 



2 Lie algebra =Sf and its 5t2-representation 

The space of all classical normal forms governed by Equation (11. 2p is denoted by This 
section provides a s[2-representation for Any column vector [/i,/2,/3]"^ is associated with 
the vector field v := /i^ + + /s^ and vice versa. Hence, v generates a system given 
by [x,y,z]'^ := [fi, f2, fs]^ ■ Besides, v acts as a differential operator on scalar functions, say 
g{x,y,z), defined by 

y(^g) ■= Ml + Ml + Ml 
dx dy dz 

Thereby, a "system" , a "vector field" and a "differential operator" are exchangeably used. For 
any vector field w, we define 

Now ^ is a Lie algebra by the Lie bracket a.dy{w) := [v, w] = vw — wv for any v,w E S£ . Any 
5I2 is a Lie algebra and is represented by a triad {M, A^, if}. In this paper they are introduced 

by 



M 



xy d xz d 



d d d 

H := 2x- y—-z—. 

ox oy oz 
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It is easy to verify [M, N] = H, [H, M] = 2M, and [H, N] = -N. The triad differential 
operators {adAf , adA?, ad//} act on the space of vector fields However, =Sf is not invariant 
under this action. Thus, any vector field from J5f generates an orbit under this action and 
the orbit is terminated when it reaches either the Zero vector field or a non-permissible vector 
field (not defined in J!f). The idea is to construct certain orbits from which we decompose the 
vector fields from ^ in terms of the orbit elements. This is done by the following definition. 
This provides a tool to use a similar approach to the method developed by Balder and Sanders 

Definition 2.1. Define 
and 

Ft ■■= ior-l<l<k, (2.1) 

Ei := ^~^l[l~^^- s.dUy' + zrEl for < / < fc, (2.2) 

e[ := ^dUy' + z^Ql iorO<l<k, (2.3) 

where ad-MV := [M, v] and ad^jf := [M, ad^^v] for any natural number n. 

These give rise to the following theorem. 
Theorem 2.2. The formulas for Fj., E\. and in terms of (x, z) -coordinates are given by 

H = + + (-1 < ^ < ^) (2.4) 

Ei = + + + (0</<^) (2.5) 

ei = x\y' + z')''-\z^-yj-), iO<l<k). (2.6) 

A straightforward calculation proves the following lemma. 

Lemma 2.3. The structure constants for the Lie algebra =Sf is governed by 

[FLf::^] = ((m+l)(A: + 2)-(/ + l)(n + 2))F^:-, 
rp/ ™i (n + 2)(m(fc + 2)-n(/ + l)) k{k + 2) 



[F^,e^] = (m(fc + 2)-n(/ + l))0'+- 

[0L,e™] = 0. 



{k + n + 2) *^+" k + n + 2 

l+m 
fc+n' 

fc+n ' 
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The algebra of first integrals for any vector field v = Uq^Fq^ + ^ ci^k-^l (when Oq ^ ^ 0) has a 
unique generator / = + ol'k^''^^ p^'^^'^^^^ in the cylindrical coordinates, see |10| Lemma 

2.1 and Proposition 2.2]. Let 

^ := spanjao ^Fq-^ + ^ a{Fl \ - 1 < / < fc, 1 < A;, 4 G m|, 

and 

T := span{eJJ + ^4et |0 < l<k,l< k,c[ G M}. 

Hence, any v & % has only phase components in cylindrical coordinates. The set of all formal 
first integrals for v G T is the algebra generated by x and y"^ + z'^, see [lOl Lemma 2.1]. Further, 
denote 

<^:={J2 ^^^^ \l + k>l,k>l,bieR]. (2.7) 

The following theorem indicates that the space of all volume-preserving and conservative 
normal forms is given by © T, see [19] for relevant results on three dimensional volume- 
preserving vector fields and their normal forms. Furthermore, it recalls that vector fields from 
(f are nonconservative. 

Theorem 2.4. The following holds. 

1. Any differential system governed by Equation (11. 2p is associated with a vector field f G =Sf 
and vice versa. 

2. The vector spaces ^ ^, and T are transversal Lie suhalgehras in S£ i.e., 

S£ = ^ ®g 

3. The Lie suhalgehra % is a Lie ideal for Thus, ^ ®^ is also a Lie subalgebra for =Sf . 

4. The algebra of first integrals for any nonzero element from % is < x-^y"^ + z^ > . 

5. Let V = Qq + aQ^^FQ^^ + J^i^i^k + '^l®!) ^ where 7^ 0. Then, there exists a unique 
formal first integral f = J^^ifi (i^odulo scalar multiplications) such that the algebra of 
first integrals for v is (/), where fl := x^^^{y'^ + z'^Y^^^^. 

6. There is no nonzero first integral for any nonzero vector field from S . Nonzero vector 
fields from S are not volume preserving. 

7. The space ^ (B is the maximal vector space of volume-preserving vector fields. The 
algebra of first integrals for any v G ^ © T zs nontrivial. 

8. The flow generated by any v & % is static in the amplitude and x- coordinates. In other 
words, it can only be dynamic in the phase coordinate. 

9. For any v & ^ ® S" , the generated flow is static in the phase coordinate. 
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Proof. The space ^ is defined such that claim [T] holds. The proofs for [2] and [3] follow Lemma 
12.31 The claims HI |5] and [6] are proved by pXil Lemma 2.1], [10, Proposition 2.2] and pT| 
Theorem 2.3], respectively. Since =Sf = © ^ © T, the proof of [7] is straightforward. Claim 
IS is true because T is generated by all vector fields from =Sf that have (only) nonzero phase 
components. Since ^ (BS' is transversal to the Lie algebra T, the flow generated by v G ^ © 
is static in phase coordinate. □ 

Theorem 2.5. The space ^ ® is a Lie algebra and there exists a Lie-isomorphism h to a 
proper Lie suhalgehra =Sfb of the Lie algebra generated by all two dimensional Bogdanov-Takens 
singularities. 

Proof. Denote for the Lie algebra generated by all vector fields of the form 

d 

w:= -x— + h.o.t.. (2.8) 
dy 

Notations are chosen such that they match with those of [5J. Define a map 

h:^®^^ ifBT, (2.9) 

governed by 

where A\ and are defined by Baider and Sanders [5l Equations 3.6a and 3.6b]. Comparing 
the two sets of structure constants, /i is a Lie-monomorphism and the claim follows by := 
h{^®<g). □ 

Remark 2.6. The claim [5] from Theorem 12.41 implies that fl{x,y,z) := x'^^(y^ + 2^)^"'+^ is 
associated with in =Sf and from [5] the monomial hl~''{x,y) := y'-+^x''-i+^ is associated with 
A^~^ in ^fc. The monomial f^'^^ = x^^"^ does not produce a permissible vector field in 
However, /i^^ = y^^"^ is associated with A'^^ . Here the s[2-orbits given by (12. ip are terminated 
when they lead to non-permissible vector fields {i.e., I is increased to k+1), while the s[2-orbits 
given by [5| Equations 3.6a and 3.6b] are terminated by zero {i.e., A'^'^ = 0). 

For bifurcation analysis of Hopf-Zero systems, there is a common practice to ignore the 
phase component. Then, study the dynamics of the two dimensional reduced system, e.g., see 
[21 [HI [ini [T71 HH]. The following theorem claims that the flow of the reduced planar system 
can be embedded into the flow generated by a Bogdanov-Takens singularity. 

Theorem 2.7. For any Hopf-Zero classical normal form, there exists a Bogdanov-Takens 
singularity such that the flow generated by the reduced planar system (Hopf-Zero) is embedded 
into the flow generated by the Bogdanov-Takens system. 

Proof. Consider the vector fields from =Sf in cylindrical coordinates. Now for any v G =Sf , let 

V := V + V, where v G T and v E Q) (S'. 
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Define w := h{v) G =2^ and i) as the two dimensional vector field obtained from v. Denote 
nonnegative real numbers by M^. The changes of variables 

/i0„ : M X M+ ^ M+ X R C M^^ 
heoix,p) = {x,y), x{x,p) := andy(a;,p) := x, (2.10) 

is a homeomorphism transforming v into w. Hence, h embeds the flow of v into the flow of w. 
This completes the proof. □ 

Notation 2.8. Throughout this paper we use Pochhammer A;-symbol notation, that is, 

fc-i 

{at:=\{{a + ih), (2.11) 

i=0 

for any natural number k and real number h. Further, for integer numbers m, n, A^, we denote 

m=^n (2.12) 

when there exists an integer k such that m — n = kN. 

Assume v = vq + vi + V2 + ■ ■ ■ , where Vk G and Jtfk denotes homogenous vector fields 
with grade k for any k. Let 

Sn [Sn,Vo]. 

Next we inductively define 

: ^„ X kerL'^-i'^-i ^ ^„ 

Af-l 

(y„,F„_i, . . . ,y„_iv+i) ^[y„_i,t;i]. 

i=0 

There exists a complement space C"'^ such that imL^ © (^n,Af _ rjj,^^ vector field w = 
Wq + Wi + W2 ■ ■ ■ is called A^-th (infinite) level normal form when Wn G C"'^ {wn G C"'"") for any 
T2,. Then, there exist invertible transformations such that v can be transformed into the n-th 
and infinite level normal form. 

3 The second level normal forms 

In this section we compute the second level normal form of the system (11. ip . 

Lemma 3.1. There exist invertible transformations transforming v^^^ into the second level 



normal form 



,(2) 



V 

k=l 

where uq = Un^ . 



e° + a,F,-' + J2{akF^ + hkE^ + c^B^), (3.1) 
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Proof. Define the grading function by S{Fl.) = = ^(0^) = k, then proof readily follows 

[Fo\ Ei] = lEi'\ [F,\ Fi] = 2{l + l)Fl~' and [F,\ &[] = 2ie\'\ □ 

The transformation Lp given by 

[x(t),p(t),0(t)] = i^{x{t),p{t),e{t)) := [x,p,e-tl (3.2) 

where G and 9 are denoted for the new and old phase variables, respectively. The changes of 
variables omits Gq from the vector field v'^'^\ Once the simplest normal form computation 
is exhausted, the map (/?~^ adds Gg back into the system, also see [HI Theorem 4.1] and [SOf 
Lemma 5.3.6]. This idea substantially reduces the computational burden and was first brought 
to our attention by James Murdock. 

After eliminating Gq, we may change ao to any arbitrary number ao via the time rescaling 
t := where t and r denote the old and new variables associated with time, respectively. 
Therefore, without the loss of generality we may assume 

oo 

v^^) = F,' + Y.{a,F^ + hEl + c,Ql). (3.3) 

fc=i 

Remark 3.2. All nonlinear Hamiltonian and Eulerian terms that may be appeared in the first 
level normal form of Bogdanov-Takens singularity in |5] are of the form A^^ and while these 
terms do not belong to =Sffe, see Theorem l2.5l Therefore, the computations in hypernormalization 
steps in this paper are essentially independent from |5] despite similarity of the procedures and 
the existence of a Lie-isomorphism between and ^ ® S . 

For hypernormalization we need to use a normal form style. A style is a rule on how to 
choose complement spaces in the normal form computation. We use the sl2-style for the second 
level normal forms. The st2-style in [5] states that the only nonlinear terms from kerad^-i can 
stay in the classical normal forms. Those are A^"^ and B\ for any /c, see [5]. Lemma [3.11 proves 
that the terms F^,E\, and G^ (see Equations fl2.ip -f l23|) ) may stay in our second level normal 
form system. These are described by ker dAp\ , if we would define F^. Indeed, we embed =Sf into 
a bigger algebra [i.e., as a proper Lie subalgebra), where the bigger algebra contains F^. As 
far as algebra is concerned, this is simply performed via a simple generalization of Lemma 12.31 
However, this extension may not have a justification in nonlinear dynamics; i.e., assigning an 
ODE system or a near-identity transformation to those extra generated vector fields (beyond 
^) is not permissible. Indeed, from the bigger algebra, F^ is only used in definitions of given 
by Equations (14. 4p . (I4.18P and fl4.30p . The map ^ is merely an instrument that substantially 
facilitates the computation. The original ideas behind definition of ^ and our gradings come 
from [5]. 

Once the second level normal form is calculated, the sl2-style is extended to a formal basis 
style. Formal basis style uses an order on, the second level, normal form terms to distinguish 
between alternative terms for elimination. Here, we give priority of elimination to conservative 
terms over nonconservative terms of the same grade in style I and vice versa in style II, see 
[T2l [T3l [211 122] for further information on formal basis style. Let 

SS^ = {Fl, E;:^,QZ\-l<l<k,0<m< n,6{Fl) = 5{E^) = 5{Q^) = N} 
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be a basis for the vector space and 5 denote for a grading function. Two formal basis styles 
are defined through the following orderings on and are used in this paper. 

• Style I: E]^ -< -< 6* for any I, k,m,n,i, j. Furthermore, -< F™, E^. -< and 

-< 0^ ii k < n. This gives a priority for elimination of i?[-terms rather than F^-terms. 

• Style II: F^ -< E[ -< Q] for any /, m, n, i,j. Besides, ^ F™, F[ ^ F™ and 6^ ^ 6™ 
if A; < n. The priority of elimination is with F^-terms rather than F^-terms. 

Throughout this paper we assume there exist ai ^ and bk 7^ for some / and k. Define 
r := min{z \ Ui ^ 0,i > 1}, s := min{j | bj 7^ 0, j > 1}, and p := min{j | cj 0,j > 1}. (3.4) 
Then, for hypernormalization we consider the following three cases: 

(i) r < s, (ii) r > s, (Hi) r = s. (3.5) 

4 The simplest normal forms 

In this section we compute the simplest normal forms for the three cases (13. 5p in the following 
three subsections. 

4.1 Case i: r < s. 

Assume 

00 00 00 

^(2) := F,-' + a.F; + J] a,F,^ + b,E', + ^ c.G^ 

fc=r+l k=s k=l 

and define 

6{Fl) = 6{Ei) = r{k-l) + k and 6{Qi) = r{k - I) + r + k + 1. (4.1) 

Note that we need to update this grading in the process of normal form computation, see 
Equation fHl5|) . 
The vector field 

:= Fo"^ + a,F; (4.2) 

plays an important role in further normalization of v^'^\ For any arbitrary a G F, a > 0, through 
changes of variables 

,asign(a.r)x 1 /Q;sign(ar). 1 ^ , ^ 

t ■= ( ^ ^ ^ )—^T, X := ( — ^ ^ )~'X, y := Y, and z := Z, 

we can change into asign(ar). Thereby, without the loss of generality we assume that in 
Equation ( 14. 2p is a non-algebraic number. When it comes to practical normal form computation 
using a computer, this assumption is not valid; computers do not recognize irrational numbers. 
However, this does not violate the practical computations, see Remak 3.2] for our reason 
in details. 
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Lemma 4.1. The (r + l)-th level normal form of U.l\) is 

v(r+^) := F,-' + arP; + hsEl + ^ a.F^ + h,El + c,Ql (4.3) 

fc>r fc>s fc>l 

where in 

• Style I, ak = for k =2(r+i) r — 1 and k =2(t-+i) r, and bk = for k =2(r+i) —I, while 
Cfc = for k =2(r+i) -1- 

• Style II, ttfc = for k =2(^+1) r — 1 and k =2(^+1) r, and for k =2(^+1) —1, Ck = for 

k =2(r+l) — 1- 

Proof. It is easy to see that v^^-* G kerad^i. Thereby, in order to compute the higher level 
normal forms we follow [4] and define 

:= ad^i oadiF^. (4.4) 

Then, 

^(i^,) = -A{l + l){k-l + 2)Pi + 2a,ik-l + l){r-l + {k-l){r + l))PlXi+\ 



, 2aJk + 2)(k-l)((k-l)(r + l) -I) 
^(i?i) = -Alik-l + l)Ei + ^ r + k + 2 

2arr{r + 2){k-l + l) 



r + k + 2 

I) = 4i^i_k~l)Q[ + 2aril-k){lir + 2)-kir + l))Q:Xi^\ 
Given the ordering for the matrix representation of ?f is lower triangular. Thus, 

ker(^) = sp3.n{p-\£l T^\keN}, (4.5) 
where P^\P^,%^ are defined in [lOl Equation 3.4]. In addition 

cO ._ ttr'^jk + 2){k)"^2 pm(r+l) \^^m^ j^Hr+l) 

m=0 ^ ^ m=0 



a^"^(/c + 2)(A;)™2 

m=0 

where S{£^) =(r+i) 0, cq = 0, and 



{r + 2{k-m){r + l))e^ r(r + 2)(A; + 1)(A;)™ 

~ 2((m+l)(r + l) + l) m!((m + l)(r + l) + l)(mr + 2(A;+l))^' 

By [TOl Lemma 3.4] we have 

[¥,..T!:] = 0, (4,7) 

|F„7?i = M!:±i)|!w^e^<-;;-, (4.8) 
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On the other hand we have 



[¥^,£l] = 2e,,^,[{2k + l){r + l) + l)F^i;^ll]l:, (4.10) 



,0 1 ^ ar''^\2k + lfT{{2k + l)(r + 1) + 2)(r + l)(2fc + 1) (2,+i)(.+ih. 

'2fc+lJ 2fc(r + 1) + 2r + 3 ^(2fc+l)(r+l)+r 



+Ae,,^,{2{k + l)(r + 1) + l)<:;)|t^;)Vt- (4-11) 

The sequence 64^+3 and 64^+2 are positive, while e^k and 64^+1 are negative. Equations (l4.7H4.9p 
imply that 

kir+l)+r 2fc(r+l)+r-l , 2fc(r+l)+r 
^fc(r+l)+r' -^2fc(r+l)+r-l ^nu -T 

are simplified in the (r + l)th level of hypernormalization. By Equations f l4.10fff?TT]) . the terms 
-^{2fc+ij(r+ij+r eliminated from the (r + l)th level normal form in style I, while in style II 
the terms F^-^l^^^^l-^^^^^^ are simplified from the system. This completes the proof. □ 

Note that the generator J^l^^^ for ker(?f ) was used to remove 

©SSlIlrt' in m while here 



2fc(r+l)+r+g 

we use it to eliminate the term i?2fc(r+ij+r+s ^2k{r+i)tr+s style I or II, respectively. 
Lemma 4.2. For any natural numbers m and n, there exist £^ such that 

n^-nT-fffji _ _|_ j2 — 2r) ((m — n)(r + 1) + m + 2)" ™ / 

rr^m IP -1 I _ ^ ^ ^ 2(r+l) j^rn-rm+n 



^(n + 2)((n-m- l)(r + 1) -m- 1) „, 

r*™ F 1 -I- F"^ = — ^ - — /-2(r+l) pj-n-rm+n 

^ '■^ " 2«-™-i(ri + 2 + r(n-m-l))(m + r + 2):';["-^ 



I n— m— 1 rprn—rm+n 



+1 

rn— rm+n ' 



where the sequence gi follows go = and 

r{r + 2)(n + 2) ((n — m — l)(r + 1) — m — l)'_2(^^;^-) 

fi'i+i = 7 \" 

2'+i(n + /r + 2)(m + r + 2)J.+i (r(/ + 1) + n + 2) (^(m + 2 + (/ + l)(r + 1)) j 

((2/ + m - n)(r + 1) + m + 2) 

--^^^ — 7 (4.12) 

2(m + 2 + (/ + l)(r + l)) ^' ^ ^ 

Proof. Let 

n-m-i n-m-1 ('|2-|-2')(fn — m — lVr + 1) — m — iV 

(Jim ._ \P ^ ^ « I \^ ^ ' AV l\ I y_-2(r+l)_ pm+/r+«+l 

" ■ a.-'2'(n + 2 + /r)(m + r + 2)^^i 



((m — n)(r + 1) + m + 2) 



'^m ' /V ' / ' /2(r+l) 

■ ^ 2'+i(m + 2)(m + r + 3)' 



«=0 



n+Z 



r+1 



Then, the proof is a straightforward calculation. □ 
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Define 

ri := r, r2 := min{i | 7^ 0, z > r}, := min{j | Cj 7^ 0, j > p} (4-13) 
where Ui and cj are associated with the (r + l)-th level normal form. 
Theorem 4.3. There exist invertible transformation sending U.l\) into 



dt 2 

d£ _ _ ar(r-+l) 

dt 



a; 



where s < r2, we have 

• /or i/ie case of style I, ak+n = for k =2r+2 —1 md k =2r+2 0, anc? (3k+s = /or 

=2r+2 -{S + 1), k =2r+2 T, and 'Jk+pi = for k =2r+2 -{pl + 1)- 

• for style II , ak+n = for k =2r+2 —1 <ind k =2r+2 and k =2r+2 s,and P^+s = for 

k =2r+2 -{S + 1), 7fc+Pi =Ofork =2r+2 "(Pl + !)• 

When s > r2 the following holds: 

• Style I: ak+n = for k =2r+2 -1 and k =2r+2 0, and j3k+s = for k =2r+2 -{s + 1), 
k =2r+2 r, and 'yk+pi = for k =2r+2 -{Pi + 1) and k =2r+2 r. 

• Style II: Ofc+n = for k =2r+2 —1 and k =2r+2 and k =2r+2 s,and j3k+s = for 
k =2t+2 -(s + 1), 7fc+p, = for k =2r+2 -{pi + 1) and k =2r+2 r. 

Besides, Equations f l4.14p represent an infinite level normal form. 

Proof. For s > r2 we change the grading function on 0-terms by 

5{Ql):=k + r2-pi-l. (4.15) 

This grading puts the term G^^ with a lower grading than F^^ and E^. Then, [TOl Equation 

3.7] imphes that J-'^^^ eliminates 02fe(r+ij+r+p| for any A; G N in the r2-th normal form level. 

Now assume that s < r2 and the grading function follows Equation f l4.ip . By Theorem 14.21 
there exist state solutions S^, 5^™, for (m, n) = (^j(r + 1) + s — 1,2k + rj + s) , such that 

^ ^o^^^-^rrir (2Mr + l) + s-r) (2(, - fc)(r + 1) + . - r);;;,, . ..^^^^^^^^^ 

2(,(r + l) + . + l)(0- + l)(r + l) + . + l)^^^7 • 



^=+1 a.2'^-^+i('=+^) (2(fc + 1) + jr + s) {2{k - j){r + 1) - s)''^ 



, 21 V j /y-v- ' ' .y- ' -/\-\- J/\- • "/-2(r+l) 2A:(r+l)+r+s 

^ 2--.(2Mr + 1) + s + 2) (0- + l)(r + 1) + sf^^^ ^^ir+.,+r.s^ 
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where 

X-fc+l _ ^ m A + A pm(r+l)-l 

m=0 ^ ' 

and (yf; follows Equation f l4.12p . Since ar is changed to be a non-algebraic number, for any G N 
either -E2fe(r+i)+r+s ^2k{r+i)+r+s eliminated from the (s + l)-th level of normalization. 

□ 



4.2 Case ii: s < r. 

For this case we assume 



OO 



V 



fc=r k=s+l k=l 

The grading function is defined by 

6{Fl) = 5{E[) = s{k-l) + k and 6{Qi) = s{k - I) + k + s + 1. (4.16) 
Denote := Fq"^ + . 

Lemma 4.4. T/ie {s + l)-t/i /ei'e/ normal form of U.l\) is 

.(^^^) = + 6.(^+^)f: + a.(^+^)F; + Y: at'' Ft + ^ h^^'^ + ^ c^M, (4.17) 

k>r k>s k>l 

where in 

• style I, b^^^ = for k =s+i 2s and k =s+i s, while cj^^^^ = for k ={s+i) s. 

• style II, a^^'' = when k =s+i 2s, b^^'^ = for k =s+i s, while cj^^^^ = if k =s+i s. 



Proof. Let 
Then, 



^ := ad(F(}) oad(E,). (4.18) 



nFl) = -4(fc - / + 2)(/ + l)Fl - + 2)JA: ^ / + 1) 



bMk-l){s + 2){k~l + l) +,+1 
^ k + s + 2 ' 

^(Fi) = -Alik~l + l)Ei-b,2ik-s)ik~l)ElX[-'\ 

'^{ei) = -Ai{k-i + i)ei + bsk{k + s)eiX'^\ 
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Then, ker(^) = span{J^^ ^^^k^Tk} where 



TT-l ._ + g + 2)(fc)^ ^rn{s+l)-l _ L m, pm(s+l)-l . 



m=l 



m=0 



- 2^ (25 + 2)-m! ^+"^^ ' ^ ^ 

m=0 ' 

for /lo = and 



+ 2)(A; - m + 2)(fc + s + 2)(A;)^-i 
~ 2™+i(s + l)'"(m)!(A; + 2 + ms)2((m + + 1) - l) 
s(m — 1) + 



2((m+l)(s + l)-l) 

Note that we have 6{T^^) =,+i -1, 5{£^) =,+i 0, d{T^) =s+i 1, and S^, = E^. On the other 
hand 

^ " ' ^ A;(2s + 2)'=(it + l)!((5 + l)(ifc+l) + l) '^"^'^'^ 

-26/+^ + 2)(. + 1) - (4-22) 

[^sA] = (4-23) 

l^-'/fcJ - 2^(s + l)^-i(A:-l)!^^+'=+^' ^ ^ 

Notice that = < 0. By an induction argument we may conclude that for any 

odd number k the number is negative and is positive for any even number k. Therefore, 
hk ^ for any k and by Equations f l4.22l - H?24l) . the proof is complete. □ 

Lemma 4.5. For any natural numbers m and n, there exist 21^ and 23™ such that 



s + 2)(n)r-iC ^sn-sm+n ^ (-6.)"-'"(n)r-(n + 2) 

sn—i 

n—m/^ ^\n—m 

n-m{^ I 1 V-rr, _L c _L n\n~m-l sn-sm+n ) 



~^[^ni^s\ — yn-n2n-m ^sn-sm+n^ 2"~™(m + 2)"~{"^"'" sn-sm+ni 



fpm _i_ roi"i m ] V" '"l—s psn—sm+n 



where — 1 < < 0, anc? 

^ (m+2):;r'-("^+^+2):;r(«+2) 



((m + s + 2)(m + 2))_ 



n— m— 1 
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Proof. The proof is straightforward by 

n ■ 2'+i(m + l)(m + s + 2)i+i '^+'^ 

"y.-^ {s + 2){ny-'{{n + ls + 2){m + 2)i^, -{n + 2){m + s + 2)1^,) 



2'+i(n + /s + 2)(m + 2)^+\ ' 

□ 

We only use style II in the following theorem. Assume that the (s + l)-th level coefficients 
Oj and bj are nonzero for some i,j > 1. Define 

Si := s, S2 := min { j I j > s} , r := min {i | 7^ 0, « > l}. 

Theorem 4.6. There exist invertible transformations sending into the {r + s + 2)-th level 
normal form 

^ = ^ _ ^^x^^/ + ife.x^^ - I Er=i K+r(A: + r + 1)0;'^ + + z ^^=1 7fc2:^ (4 

^ = -y- ^^x^z + \bsX^z - f {(^k+rik + r + l)x^ + Pk+sX')x>' - vY^Zi lkx\ 

where a/c = for k =s+i 2s, /3fc = for k =s+i s, 7^ = for k =s+i s, and 
if S2 < r, we have (3k = for k = s"^ + S2 + s, 
when S2 = r, = for k = s"^ + S2 + s, 
for S2 > r, al^~^^^ = for k = s"^ + r + s. 

Furthermore, Equations fl4.25p is indeed the infinite level normal form. 
Proof. For 82 < r, by Lemma 1^31 there exists a ^^+53 ^^^'^ ^^^^ 

rais2 Iff 1 I rz?o 771521 _ (~^g)^(-^2 ~ ■s)(-S2)s ^52^,^^^^ 

So, in both styles we can remove -£"^2+^2+5 f^o^ the system. For S2 = r, by Lemma 14.51 there 
exist 21^^52 si^ch that 



\{s2 + l){s2 + s + 2r-\ + (., + . + 2)^(., + . + 2):;} + + 2)(.,),,,, J jE^ 

(g2 + 2)(g2):+^fls2 z^.2+S2+S 
(s2 + 2):+} ^'+^^+^' 



'+S2+S 
+S2+S 



M. Gazor and F. Mokhtari 



Hop f- Zero Singularity 17 



Therefore, F^2^g^^g can be eliminated in style II. 

When S2 > r, by Lemma H75] there exists a 21^ such that 



2(r + s + 2) V ^ (r + s + 2 

(_f,^)V(r + 2)(r + s)f 



2(r + 2):+i 

Since both coefficients are nonzero, the proof is complete. □ 



4.3 Case iii: r = s. 

Let 

oo oo 

v(^) := F,-' + asF: + h.El + ^ (a,F,^ + hEl) + ^ c.G^ (4.26) 

k=s+l k=l 

and assume that is a non-algebraic number. This section is similar to the case III of Baider 
and Sanders [1] which is the most difficult case of the three. We essentially use the idea of 
Wang et al. [151 El]; they assume a certain ratio be non-algebraic. This prevents that the 
ratio be a root for any polynomial that appears in the computations and the results readily 
follows. When the corresponding fraction is an algebraic number, both problems (the case in 
this subsection and the third case of [4]) remain unsolved. Denote 

Xs.= F^^ + asF! + hsEl, (4.27) 

where 5{^s) = and define a grading function by 

5{Fl) := 6{Ei) := s{k - I) + k and 5{Q{) := s{k - I + I) + k + I. (4.28) 

Theorem 4.7. Let ^ be a non- algebraic number. Then, there exist invertible transformations 
that sends v^^^ given by equation into the infinite level normal form system 

I = (2/2 + z') + (b, + a,)x^+i + X EZs+iiPk + afc)x^ 

^I = z+{bs- as{s + 1))^ + f EZs+i {^>^ - Mk + l))x'' + z Er=i 7.a;^ (4.29) 
f = -y + (fe. - asis + 1))^ + f EZs^iifik - au{k + l))x^ - y Y.Zi 

where in 

• style J, = for k =s+i —1 and k =^+^ 2s, and for k =s+i —1 we have 7^ = 0. 

• style J/, afc = for k =s+i —1 and k =s+i 2s, while 7fe = for k =s+i — 1- 
Proof. Define 

$f := ad(F(}) o ad(X,) (4.30) 
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and 



-2(fc-0(- 



Hence, 



where 



7^ 



I , 4a, 5(5 + 2)(A;-/ + l) 
k + s + 2 
'a,{k + 2){l-{k-l){s + l)) 



k + S + 2 



-b,ik-s)]EiX^;\ 



-2l{k-l + l)Qi+ (bsk{k + s)-2ar{k-l){l-{k-l){s + l))^e[X'+\ 



ker(^) = span{ 7^, £° | A; e N}, 



(4.31) 



fc+i 



k+l 



m=l m=l 
k k 

m=0 

E 



-1/, m pm(s+l)-l 



Of,™ jj,m(s+l) 



m=0 



m=0 



m!(2(s + l))"»-i 



5(7-,-^) ^,+1 (-1), 
5{£l) =s+i 0, 
5(7;?) ^,+1 (-1), 



and ^ = — 0- The sequences ul^ and for i — 0,-1 and m > follow the recurrence 
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relations 



-1 ass{s + 2) _^ 

K{m + l){s + l){m{s + l) + k + 2f'^ ^ 

{s + l)(2(m - 1) - k){m{s + 1) + A; + 2) + hs{ms + k){ms + k + 2) , 

^— 11 

6,(m + l)(s + 1) (m(s + 1) + A; + 2) 

- s{s + 2){k-m + 2) 

■ [m{s + l) + s){s{m+l) + k + 2) ^ ^ ^ 

ttsims + k + 2){{k- 2m + l)(s + 1) + l) - bs{s{m- 1) + A;) (m(s + 1) + A; + 2) _^ 
^ (m(s + 1) + s) {s{m + l) + k + 2) ^"^ ' 

■= (4-34) 
&4s(m + 1) + A; + 2) 



(s(m + 1) + A: + 2) ((A; - 2m)(s + 1) + s + 2) - 6,(ms + k){ms + k + 2) 



bs{s{m + 1) + k + 2)' 





m ' 



- f!^^?^""':^' ,.,„. < (4.35) 



+ l)(s + l)(s(m + l) + A; + 2) 



as{s + l){k - 2m)(ms + A; + 2) - 6Js(m - 1) + A;) (s(m + 1) + A; + 2) „ 

H — -w 

bs{m+l){s + l){s{m + l) + k + 2) 

On the other hand 

[x„ j--^] = 26/'+2(^ + 2){u-i,{s + + <2((^ + 2)(3 + 1) - i)Ei:xixi)., 
[x„£:°] = 26/'+^(t.°_,,((fc + i){s + 1) + i)F,%tt^ + + i){s + i)i?Lt,^::), 

^ " - 2'^(A;-l)!(s + l)^'-i ^ks+k+s- 

Using the recurrence relations (14 .3 2114^351) . one can compute w^^2' ""fc+i' '^l+i 

terms of a polynomial of ^ with rational coefficients. Since ^ is non-algebraic, all these four 
numbers are not zero and the proof is complete by the rules established in styles I and II. □ 



5 Examples 

In this section we provide the formulas for the first few coefficients (up to degree 3) of the 
simplest normal forms in terms of the coefficients of the original system (before classical normal 
forms). These are very useful for practical applications. Then, there are applied on the Rossler 
and Kuramoto-Sivashinsky equations to demonstrate the applicability of our results. The 
following proposition presents the formulas for the first few terms of classical normal form in 
terms of conservative-nonconservative decomposition introduced in this paper. 
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z'y^x'', (5.1) 

where Co,2,o + C2,o,o 7^ 0. Denote its cubic-truncated classical normal form by 

^« 00 + a,'F^' + aSF° + b\E[ + a\Fl + b\E\ + a^F| + blEj + c}0} + c?0? + c^0^. 

= ^(C0,2,0 + C2,0,o), (5.2) 

g(-2co,o,2 + ai,o,i + 

fi2 ■ — o (~2Co.o,3 ~ 4ao,0.2Co,l,l + 46o,0,2Cl,o,l + <3^1,0,2 + ^0,1,2 + Oo,0,2«l,l,0 + 2ao,o,2&0,2,0 

~2a2,o,oVo,2 ~ ^0,0,2^1,1,0)) 
'■— ^(— 8co,2,i — 4ci_o,2 — 5ao,i,iCo,2,o + 5ao,i,iC2,o,o ~ 5ai_o,iCi,i,o + 201^1^0^2,0,0 

— 4a2,o,o&2,o,o + 56o,i,iCi,i,o — 4&o,2,o&2,o,o — 96i,o,iCo,2,o + 56i,o,iC2,o,o 
—26i i 062,0,0 + 2ao,2,o + 603^0^0 + ^2,0 + 262,1,0 ~ 8«'0,2,oCo,i,i 
— 8a2,o,oCo,i,i + 86o,2,oCi,o,i + 862,o,oCi,o,i), 

— ^(co,o,2 + ai,o,i + 60,1,1), 

62 '■— g(— 2co,o,3 ~ 12ao,o,2Co,i,i + 126o,o,2Ci,o,i + 5ai^o,2 + 56o,i,2 + 5ao,o,2ai,i,o + 1000,0,260,2,0 
— 1002,0,960,0,2 — 560,0,261,1,0), (5-3) 

^1 •= 7^(~3co,2.i - 8ci,o,2 - 13oo,i,iCo,2,o + 13ao,i.iC2.o,o - 13ai,o,iCi,i,o + 10oi,i,oa2,o,o 
—2002,0,062,0,0 + 136o,i,iCi,i,o — 2060,2,062,0,0 ~ 216i,o,iCo,2,o + 136i,o,iC2,o,o 
-1061,1,062,0,0 + 10oo,2,o + 30o3,o,o + 1060,2,0 + 1062,1,0 - 16oo,2,oCo,i,i - 1602,0,OCO,1,1 
+166o,2,oCi,o,i + 1662,o,oci,o,i, 
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c\ := 7(-ao,i,i + bi,o,i)> 



and 

4 

_ 1. 3 1 3 1 5 5 1 

Ci :— 2 ^ ~ g^o,3,o — g'^2,1,0 + gt^3,o,o ~ Y^Oo,i,iCi^i_o + ^fli,i,oOo,2,o ^ ^2^°'^'°^^'°'° ~'~ 24*^°'^'° 

1 1 1. 2 5 2^ 1 1 2^ 5 

1, 1, "^1,2^1,7 '^/ii,2^^2\ 

~ 77700,1, l'^0,2,0 + 77700,1, l'^2,0,0 ^ 7700,2,0 ~ 7700,2,002,0,0 ~ 777Ol,0,lCl, 1,0 — -40l, 1,0 —7702,0,0 )) 

lb lb b 1/ lb / 12 

C2 7(~2ao,l,2 — 261,0,2 — 400,0,20-0,2.0 + 2(30,0,2^1,1.0 — 4ao.o,2Cl,0,l ~ '^0,1.1^ ~ '^0,1, 1^1,0.1 — 77'3l,0,l^ 

+ai,o,iVi,i + 2ai,i,o6o,o,2 - 460,0,2^2,0,0 - 46o,o,2Co,i,i - (5-4) 

Proof. The formulas are derived by explicit calculation. They are implemented by treating the 
coefficients as symbols in Maple and verified by the program. □ 

We are interested to derive formulas for the coefficients of the simplest normal forms in 
terms of the system fl5.1l) . Given the formulas fl5.2H5.4l) . we may present the formulas in terms 
of the coefficients of classical normal forms. Instead, we define the numbers ao, fli, 02, 03, 61, 62 
and Ci in terms of the coefficients of classical normal forms by 



Co 


:= ^(co,: 


ai 


:= \{c., 


hi 


■■= \{c,, 


Cl 


■■= 4(^0,: 




'^0 (. 


02 




4 



1 



Then, the simplest normal form coefficients are given in terms of Oj, 6j and Ci. 

Theorem 5.2. Assume that ai^ + 61^ 7^ 0. Then, the three-degree truncated simplest normal 
form of Equation / (5. 1\) is given by 

dec 

— = 2aop^ + (ai + bi)x'^ + (02 + 62)x^ (5.5) 
^ = ^(61 -ai)xp + i(6i -3ai)xV, 

de 
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Proof. The proof follows a symbolic implementation of the results in Maple. □ 
Example 5.3. Consider the Rossler equation 

X = -y- z, 

y = X + ay, (5.6) 

z = bx — cz + xz. 

This system has Hopf-Zero singularity for two sets of parameters. One is for parameter values 

a = c, 6=1, < < 2, (5.7) 

while the system associated with the other set of parameter values has a simple dynamics, see 
[1]. Hence, we denote the vector field associated with (15. 6p by Va- Then, Proposition 15.11 implies 
that the classical normal form of Equations (15. 6p is given by 

= e° " I ^(4«' + 7) ^0 a{2 + a^) 

° y2^(-2 + a2) ° 8(2-a2)f(-2 + a2) ' 6v/2^(-2 + a^) ' 
3(2a2 + l)a 3a(4a2 + 5) a(-l + a)(a + l) 



K2-a2)i(-2 + a2) 8(2 - a^)! (-2 + a2) ?,y/2^{-2 + a"^ 
(22a2 + 7)a ^ + 1) , (-21a^ + 12a^ - 10) 



3(2-a2)f(-2 + a2) ' 4(-2 + a^) ' 48(-2 + a2)3 

fa^ + 62a^ - 760^ - 8 



64(-2 + a2)= 



-6 



2- 



Then, 



R'r' = e'o-Fo'-F^-lF^-lEi + Ql 

and for a 7^ ±1 we have 

" " °+(2-a2)i ° +v/2^(-2 + a2)^^+ 3(-2 + a2)3 

V2 - a2(16 - 32a^ + 56a2) + 16(2 - a^)! (^2 + 1) + 3a{3a^ - 2a^ - 6a^ + 2a - 5) 



-32a-2(2 - a2 



11 
2 



V2^^(96a^ - 48 - 168a2) + 16(2 - a2)2(a2 + 1) + 3a{3a^ - 2a^ - 6a2 + 2a - 5) ^2 

n ^2 

32a-2(2-a2)- 

a(a2 + 1) 



2V2^(-2 + a2)2 



el. 
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Now we apply our results on the well-known Kuramoto-Sivashinsky equation. 
Example 5.4. Consider the Kuramoto-Sivashinsky equation 

d 

ht + -^Gih, K,K,) = 0, (5.8) 

where 

XXX ■ 

A common technique in transforming a partial differential equation into a system of ordinary 
differential equations is to assume that the solutions are as traveling wave solutions given by 
h{x, t) := u{x + fit). This gives rise to 

X = y, 
y = X, 

z = jjix — 2x'^ — y, 

for some parameter /i. For /i = 0, the system has a Hopf-Zero singularity at the origin, see 
Chang [7J. Our formulas for classical normal form of up to degree 3 is given by 

^ 96 ^ ^ 96 ^ 24 ^ 4 2 
We have Oi 7^ and 61 = 0. Thus, the infinite level normal form is 

X = + 16x^, 

p = —16xp — 31x^p, (5.9) 
e = 1. 
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